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THE FADE INTERPOLATION METHOD APPLIED TO ^-PAINLEVE EQUATIONS 

H (DIFFERENTIAL GRID VERSION) 


HIDEHITO NAGAO 

Abstract. Recently we studied Fade interpolation problems of ^-grid, related to g'-Painleve 
equations of type ' and (A 2 + By solving those problems, we could 

derive evolution equations, scalar Lax pairs and determinant formulae of special solutions for 
the corresponding ^-Painleve equations. It is natural that the g'-Painleve equations were derived 
by the interpolation method of q'-grid, but it may be interesting in terms of differential grid that 
the Fade interpolation method of differential grid (i.e. Fade approximation method) has been 
applied to the ^-Painleve equation of type by Y. Ikawa. In this paper we continue the above 
study and apply the Fade approximation method to the g'-Painleve equations of type Eg*', 

A^' and (A 2 + Ai)^'\ Moreover determinant formulae of the special solutions for q'-Painleve 
equation of type Eg*' are given in terms of the terminating ^-Appell Lauricella function. 


1. Introduction 

In this paper we eontinue [fT3]l and apply the Fade approximation method to the qr-Painleve 
equations of type Eg \ ^ and (A 2 + 

1.1. The background of discrete Painleve equations. 

In Sakai’s theory ll^ the diserete Painleve equations were elassified on the basis of rational 
surfaces connected to extended affine Weyl groups. There exist three types of diserete Painleve 
equations in the elassification: elliptie difference {e-), multiplicative differenee {q-) and additive 
difference (J-). The diserete Painleve equations of ^-differenee type are elassified as follows: 
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Here A ^ B means that B is obtained from A by degeneration. 

1.2. The background of the Fade method. 

Fade approximation/interpolation are closely related to Painleve/Gamier equations. The Fade 
method is a method for giving Painleve equations, sealar Lax pairs and determinant formulae 
of special solutions simultaneously, by starting from suitable problems of Fade approximation 
(of differential grid)/interpolation (of differenee grid). In 1(311 Y.Yamada has applied the Fade 
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method to continuous Painleve equations of type Pyi, Py, Pw and Gamier system by using 
differential grid (i.e. Fade approximation). 

The Fade method for discrete Fainleve equations has been applied to the following types: 
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It is natural that the continuous/discrete Fainleve equations were derived by the interpolation 
of differential/difference grid respectively. Here it may be interesting to note that the Fade 
approximation method of differential grid (i.e. Fade approximation) has been applied to the 
type in BHI. In this paper differential grid is applied to type ^ 

q-{A2+A,i^\ 

Remark 1. On the key points of the Fade method 

There are two key points to apply the Fade approximation/interpolation method llH |T3l |T7l 
|33l . The first key point is the appropriate choice of approximated/interpolated functions (see Ta¬ 
ble (12.21) and Remark H]). The second key point is to consider two linear ^-difference three term 
relations (12.61) satisfied by the error terms of the Fade approximation/interpolation problems. 
Then the error terms can be expressed in terms of special solutions of ^-Fainleve equations. 
Therefore the ^-difference relations are the main subject in our study, and they naturally give 
the evolution equations, the Lax pairs and the special solutions for the corresponding (^r-Fainleve 
equations. □ 

Remark 2. On a connection between the Fade method and the theory of semiclassical 
orthogonal polynomials 

The connection between semiclassical orthogonal polynomials (classical orthogonal poly¬ 
nomials related to a suitable weight function) and Fainleve/Garnier systems has been demon¬ 
strated in IfTTII . It has been shown that coefficients of three term recurrence relations, satisfied 
by several semiclassical orthogonal polynomials, can be expressed in terms of solutions of 
Fainleve/Garnier systems (see [|Il[l5l|T8l|2S|271|^|23 for example). Thus there exists a close 
connection between the Fade method and the theory of semiclassical orthogonal polynomials. 
Namely, using both approaches, we can obtain the evolution equations, the Lax pairs and the 
special solutions for the corresponding Fainleve/Garnier systems. (The theory of semiclassical 
orthogonal polynomials is more general and the Fade method is simpler. For example their 
relation was briefly proved in OOl .) □ 

1.3. The purpose and the organization of this paper. 

The purpose of this paper is to apply the Fade approximation method to type \ ^ 

A 4 ^ and q-{A 2 -l-As the main results given in Section [3l the following items are presented 
for each type. 

(a) Setting of the Fade approximation problem, 

(h) Contiguity relations. 
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(c) The Painleve equation, 

(d) The Lax pair, 

(e) Speeial solutions. 

This paper is organized as follows: In Section!^ we explain the Fade approximation method 
applied to the ^-Painleve equations, namely the methods for the items (a)-(e) above. In Seetion 
|3]we present these main results for type q-E^^\ ^ q-{A 2 + In Seetion|4]we 

give a summary and diseuss some future problems. 


2. Fade approximation method of differential grid 

In this seetion we explain the methods for deriving the items (a)-(e) in the main results given 
in Seetion |3l These eontents of the items (b)-(d) below (i.e. Subseetion 12.2[ 12.31 and 12.41) are 
almost the same as the items (b)-(d) in Seetion 2 of iTTSl . 

2 .1. (a) Setting of the Fade approximation problem. 

Let us eonsider the following approximation problem (of differential grid): 

For a given funetion Y{x), we look for funetions Pm{x) and Q„{x) whieh are polynomials of 
degree m and n e Z>o, satisfying the approximation eondition 

(2.1) Y{x) = ^^^ (mod 

Qn{x) 

We eall this problem the ’’Fade approximation problem (of differential grid)”. Then the funetion 
Y{x) is ealled the ’’generating function” (because Y{x) generates the coefficients in power 
series (12.131) in the item (e) below, i.e. Sub section [23]), and the polynomials Pmi^) and Qn{x) are 
called ’’approximating polynomials” respectively. The explicit expressions of the polynomials 
Pm{x) and Q„{x), which are used in the computations for the item (e) above, are given in the 
formulae (12.141) and (12.201) (see the item (e) below). 

Remark 3. On the common normalization factor of the polynomials F„,(a) and 2 „(a) 

The common normalization factor of the approximating polynomials Pmix) and Qn{x) is not 
determined by the condition (12.11) . However this normalization factor is not essential for our 
arguments, i.e. the main results in Section [3] (see Remark[5|). □ 

Fix a complex parameter (0 < |^| < 1). Let a,-, h, e be complex parameters. We establish 
the approximation problems (12.11) by specifying the generating functions Y{x) as follows: 
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a\a2a^q o) 

Here - = 1 is a constraint for the parameters in the case q-E). , and the (^-shifted 

factorials are defined by 


(2.3) 


./-i 

(ai,a2, ■ • ■ ,ai-,q)j := ]~^(1 - aiq’"){\ - a2q'') ■ • ■(! - atq'') 


k=Q 


and the ^-HGFs (i.e. the (^-hypergeometric functions [0) defined by 

X (au---,ak;q), r. . . mi!+/-<= 

(2.4) k‘Pi 


"" :=y k-D-#] 

bu ■■■, bi I ^ (bi,--- ,bi,q;q),^ J 




with ( 2 ) = 5(5 - l)/2. 

The hypergeometric solutions to the ^-Painleve equations of type ^ and (A 2 + 

Ai)'^^^ were given in terms of the (?-hypergeometric functions 3 ^ 2 , 2 ^ 1 , 2 < 4 >i, and ii/?i respectively 
in [[^. The functions Y{x) in Table (12.21) generate the coefficients pk in power series (12.131) in 
terms of the terminating ^-hypergeometric functions 31/72 (13.131) (see Remark|7]), 2 <P\ (13.251) . 2 <P\ 
(13.361) . and lipi (13.471) . 


Remark 4. On the choice of the generating functions Y{x) 

One may wonder how the generating functions Y{x) are appropriately chosen. However there is 
no theoretical choice of the functions Y{x) in the Fade approximation method as far as we know. 
We chose the functions Y{x) to type E^^\ A^ ^ and (A 2 + Ai)^^^ by a extension and reductions 
from the generating function to type given in HI. □ 

Let us consider yet another Fade problem where some parameters a,, bu m and n in the gen¬ 
erating functions Y{x) are shifted. The parameter shift operators T are given as follows: 
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Here the operators T are called the ’’time evolutions”, because they specify the directions of the 
time evolutions for ( 2 '-Fainleve equations. 


2.2. (h) Contiguity relations. 


Let us consider two linear three term relations: L 2 {x) = 0 between y(x),y(( 2 'x),y( 3 c) and L 3 (x) = 0 
between y{x), y{x), y(xlq) satisfied by fundamental solutions y(x) = Pm(x), Y{x)Q„{x), where L 2 
and L 3 are given as expressions 
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Then the linear relations L 2 = 0 and L 3 = 0 are called the ’’contiguity relations”, and the conti¬ 
guity relations are the main subject in our study. Here for any object F the corresponding shifts 
are denoted as F := T{F) and F := T~^{F), and the shift operator T acts on parameters given 
in Table (12.51) . 


We show the method of computation of the contiguity relations L 2 = 0 and L 3 = 0. 


SetyCx) : = 


Y{x)QAx) 


and define Casorati determinants Di{x) by 


(2.7) Di{x) := det[y(ji;),y(< 5 rA)], D 2 {x) := det[y(T), y(ji;)], D 2 ,{x) := det[y(< 5 rT), y(A)]. 
Then the expressions (12.61) can be rewritten as follows: 

L 2 {x) oc Di(x)y(x) - D 2 {x)y{qx) + D^(x)y(x), 

L^ix) oc Diixlq)yix) + D 2 {xlq)y{x) - D 2 (x)y(x/q). 

Define basic quantities G(x), K{x) and H(x) (e.g. (13.51) . (13.171) 1 by 

(2.9) G{x) := Y(,qx)/Y(x), K{x) := Y(x)/Y(x), H{x) := L.C.M(Gde„(x), Fden(T)). 


Here Gden(T) and Gnum(.T) are defined as the polynomials of the denominator and the numerator 
in G{x) respectively, and K^^nix) and K^amix) are similarly defined. For example in case of q- 

Fg , Gden('*l) ~ Oi=l(l GiX\ Gnum(T) — ]~[j_3(l biX\ Fden(x) — 1 CliX, Fnum(T) — 1 biX 

(see eq. (13.51) 1. Substituting these quantities into the determinants (12.71) . we obtain the following 
expressions: 


D,{x) = 
(2.10) F>2 (x) = 


Y(x) 

(-^den(T) 

Y{x) 

^den(^) 


{G]^am(.X^Pm(.X^Qn(.^X^ (-^den(^)^m(9'^)2«(^)} ’ 

^num(.X^Pm(^x}Qi^(^x} Fden(T)F,j,(x)2H(T)| , 


D3{x) 


Y{X) 

H(x) 


Hjx) 

Fden(T) 


Fnum (^x)Pffi {fl X)Qfi{^x) 


H{x) 

6^den(T) 


f^num {x)Pm{x)Qn{qx) 


Using the approximation condition (12.11) and the form of the basic quantities G(x), K{x), H{x) 
(e.g. eqs. (l3.5l) . (13.171) 1. we can investigate positions of zeros (e.g. x = 0) and degrees of the 
polynomials (e.g. Gnum(x)Pm(x)Qn(qx) - Gden(x)P,„(qx)Qn(x)) within braces { } of the expres¬ 
sions (12.101) . Then we can simply compute the determinants D,(x) (e.g. eqs. (l3.6l) . (13.181) 1 except 
for some factors such as 1 - fx, 1 - xjg and c, in D,(x), where /, g and c, are constants with 
respect to x (see Remark!^. In this way we obtain the contiguity relations L 2 = 0 and F 3 = 0 
(e.g. eqs. dXTl) . (13.191) 1. 


Remark 5. On the gauge invariance of the product CqCi 

When the common normalization factor of the approximating polynomials Pm{x) and Qn{x) is 
changed, an x-independent gauge transformation of y{x) is induced in the contiguity relations 
L 2 = 0 and L 3 = 0. Under the x-independent gauge transformation of y(x): y(x) Gy(x), 
the coefficients ofy(x),y{x/q),y(x) and y(x),y(x),y(x/q) in L 2 = 0 and F 3 = 0 (12.81) change as 
follows: 


( 2 . 11 ) 


(Di(x) : D2(x) : D,(x)) ^ {GD,{x)JG : D2(x) : D^ix)) 

(Di(x/^) : D^ixlq) : D2(x)) 1 -^ {GDi{xlq)IG : D^ixlq) : D2{x)). 
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The coefficients Co and Ci in L 2 = 0 and L 3 = 0 (e.g. eqs. (l3.7L (13.191) t are defined as the 
normalization factors of the coefficients of y(x) and y(x) respectively. Then Cq and Ci change 
under the gauge transformation, but the product CqCi is a gauge invariant quantity. Moreover 
Co and Ci do not appear in the final form of the (?-Painleve equations (e.g. eqs. (l3.8L (13.201) 1. □ 

Remark 6. On two meanings of variables /, g and parameters m, n 

We use / and g with two different meanings. The first meaning is constants (i.e. special so¬ 
lutions) / and g which are explicitly determined in terms of parameters a, , bi, m and n by the 
Fade approximation problem (e.g. eqs. (l3.6l) . (13.71) . (13.121) . (13.181) . (13.191) . (13.241) 1. The second 
meaning is generic variables (i.e. generic solutions) / and g apart from the Fade approximation 
problem (e.g. eqs. (l3.8l) . (13.111) . (I3.20L (13.231) 1. namely / and g are unknown functions in the 
^-Fainleve equation. In the items (c), (d) (resp. in the items (b), (e)) we consider / and g in the 
second meaning (resp. in the first meaning). 

Similarly we use m and n with two meanings. In the first meaning m and n 6 Z>o are integer 
parameters (e.g. eqs. (l3.3l) . (13.61) . (13.71) . (13.151) . (I3.18L (13.191) 1. In the second meaning m and n 
e are generic complex parameters, namely q'" and are replaced by generic parameters ao 
and bo respectively (e.g. eqs. (l3.8l) . (13.111) . (I3.20L (13.231) 1. In the items (c), (d) (resp. in the items 
(a), (b), (e)) we consider / and g in the second meaning (resp. in the first meaning). Then the 
result of the compatibility of the contiguity relations L 2 = 0 and L 3 = 0 also holds with respect 
to the second meaning. □ 

2.3. (c) The (^-Painleve equation. 

Let us consider generic variables /, g and generic parameter ao, bo as in the second meaning in 
Remark 0 we can derive the ^-Fainleve equation as the necessary condition for the compat¬ 
ibility of the contiguity relations L 2 = 0 and L 3 = 0 (e.g. eqs. (l3.7L (13.191) ). Computing the 
compatibility condition, we determine three variables g, f and CqCi. Expressions for variables 

g and / are obtained in terms of variables / and g. An expression for the product CqCi is 

obtained in terms of variables /, g and / (and hence in terms of variables / and g). 

The first and the second expressions are the ( 2 '-Fainleve equation (e.g. eqs. (l3.8L (13.201) 1. The 
third expression is a constraint for the product CqCi (e.g. eqs. (l3.9l) . (13.211) 1. 

2.4. (d) The Lax pair. 

Let us consider two linear three term equations for the unknown function y(x): Li{x) = 0 
between y(qx),y(x),y(xlq) and L 2 (x) = 0 between y(x),y(qx),y(x), where Li and L 2 are given 
as expressions 

Li(x) = Ai(x)y(xlq) -H A 2 (x)y{x) + A 2 (x)y(qx), 

L 2 {x) = A 4 (x)y(x) + A 5 {x)y(x) + Ao(x)y(qx). 

The linear three term equations L\ = 0 and L 2 = 0 (12.121) are called the ’’scalar Lax pair”, when 
the compatibility condition of the linear equations Li = 0 and L 2 = 0 ( 12 . 121 ) is equivalent to a 
(^r-Fainleve equation. 
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We present how to eompute the sealar Lax pair. Similarly to the item (e), let us eonsider 
generie variables f,g and generic parameter aQ,bo as in the second meaning in Remark 0 The 
Lax pair L\ = 0 and L 2 = 0, which satisfies the compatibility condition, is derived using the 
results of the items (a)-(c) as follows: The Lax equation L 2 = 0 (e.g. eqs. (l3.11l) . (13.231) 1 in the 
item (d) is the same as the contiguity relation L 2 = 0 (e.g. eqs. (l3.7L (13.191) 1 in the item (b) 
under an x-independent gauge transform of y{x) and changes of parameters. We can obtain the 
Lax equation L\ = 0 as follows: First combining the contiguity relations L 2 = 0 and L 3 = 0 
(e.g. eqs. (13.7l) . (13.191) 1 under generic variables f,g and generic parameters ao,bo, one obtains 
a linear equation between the three terms y(qx),y(x) and yixjq) (See the figure below), whose 
coefficient functions depend on the variables f,g,f,Co and Ci. However the variables Cq and 
Cl appear through the product CqCi. Therefore expressing / (e.g. eqs. (l3.8L (13.201) 1 and CqCi 
( e.g. eqs. (13.91) . (13.211) 1 in terms of / and g only, one obtains the Lax equation Li = 0 (e.g. 
eqs. dLTO . (IL23]) 1. 


y(x/q) y(x) 



y(x/q) Li{x) y(x) y(qx) 

Then, under generic variables f,g and generic parameters ao,bo, the ( 2 '-Painleve equation 
(e.g. eqs. (l3.8L (13.201) 1 is sufhcient for the compatibility of the Lax pair Li = 0 and L 2 = 0 (e.g. 
eqs. (l3.11L (13.231) 1. The similar proofs have already been given in lfgl [3n[32l . 

2.5. (e) Special solutions. 

By construction, expressions for constants / and g as in the first meaning in Remark!^ give a 
special solution for the i^r-Painleve equation. We present how to compute determinant formulae 
of the special solutions / and g. 

We derive the formulae (12.141) which are convenient for computing the special solutions / and 
g. We can assume the power series 

00 

(2.13) Y{x) = J]ptx\ po = l, Pi = 0 (/<0) 

k=0 

since the generating functions Y{x) in Table (12.21) are holomorphic near x = 0. Then for each 
type of the functions Y(x), the coefficients pk (e.g. eqs. (13.251) . (13.361) 1 are given respectively as 
the terminating cases of ^-hypergeometric functions defined by eq. (12.41) . 

For a given function Y{x), the polynomials P„j(x) and 2„(x) of degree m and n for the ap¬ 
proximation condition (12.11) are given by the following determinant expressions: 

m n 

(2.14) Pffi(^X^ ^ 5 QnC^) ^ 'j •^((m+iy,m"^')( 

i=0 i=0 

where sa is the Schur function defined by the Jacobi-Trudi formula 

(2.15) := det(p2._i+^')!j=i. 
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We show the derivation of the expressions (12.141) as follows: The approximating polyno¬ 
mial Qn{x) satisfying the eondition (12.11) ean be given as the seeond expression of eg. (12.141) by 
Cramer’s rule. The approximating polynomial Pm{x) satisfying the eondition (12.11) is given by 
the following eomputation: By using the relation 


(2.16) v:"T(x) = 

k=Q k=Q 


we have 
(2.17) 



Pm 

Pm-Yl 

Pm+n 

CO 

Pm 

Pm-\-{ 

Pm+n 

Y(x)Qn(x) = 

Pm—n-Y\ 

x^Y{x) 

Pm 

• ■ • xY{x) 

Pm+l 

Y{x) 

= Z 

k=0 

Pm~n+l 

Pk~n 

Pm 

• • • Pk-\ 

Pm+l 

Pk 


Here we note that 


m+n 

(2.18) ^ S(m\k)^ = 0. 

A:=m+1 

Substituting the relation (12.181) into the expression (I2.17L we obtain 

m oo 

(2.19) Y(x)Qn(x) = (E- Z 

k=0 k=m+n+l 


oo 

^ ,k)^■ 

k=0 


Henee the desired polynomial Pmix) is given as the first expression of the formulae (12.141) . 

Furthermore the polynomials P,n{x) and Qnix) in the formulae (12.141) ean be expressed in 
terms of a single determinant as 


( 2 . 20 ) 


P,„{x) = ^ (-x)”5((„+i)»)| 


\Pi^Pi-x V/-r 


The formulae (12.141) and (12.201) have already appeared in 0^ . 

Then we apply the general results deseribed above to the ease N = I (13.151) and N = 2 (13.31) 
of the funetion Mx) := OtV , whieh ean be written as 

(biX)oc 


( 2 . 21 ) 


il/(x) = exp 


/ K-ai 




^) = Z 


PkX^. 


k=Q 


We note that this kind of expression (12.211) has already appeared in 


We show the method of eomputation of the speeial solutions / and g. 

The expressions for the speeial solutions / and g ean be derived by eomparing the determinants 
Di(x) in eg. (12.101) and Di{x) (e.g. egs. (l3.6l) . (13.181) 1 in the item (b) as the identity with respeet to 
the variable x and applying the formulae (12.141) and (12.201) . 

For example the eomputation for the ease is as follows: Substituting x = 1 jui (i = 1,2) 
into the determinants T>i(x) in eg. (12.101) and Di(x) in (13.61) respeetively, we obtain an expression 
for the speeial solution / in the first eguation of eg. (13.121) by eomparing the two expressions for 
Di(x) and applying the formulae (12.201) . Similarly substituting x = 1/h, (i = 2,3) into the 
determinants D^(x) in eg. (12.101) and D^(x) in eg. (13.61) respeetively, we obtain an expression for 
the speeial solution g in the seeond eguation of eg. (13.121) by eomparing the two expressions for 
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D^{x) and applying the formulae (12.201) . 


3. Main results 


In this seetion for eaeh ease q-E^^\ q-^ 4 ^ and q-iA 2 + we present the results 

obtained through the method, whieh was explained in Seetion 
We use the following notations: 


aia2 ■ • ■ anlb\b2 • • - bn 


Cl\Cl2 ' ’ ’ Cln 

b\b2 ■■■bn 


(3.1) 




Tai(F) ■= F\a,-.qa,, != F\a^^a./q 

for any quantity (or funetion) F depending on variables a, and bi, and by definition (12.151) . the 
Sehur funetion is expressed as 

Pm+n—i 


Pm 

where the element pk is defined in the power series (12.131) . 


(3.2) 




Pm Pm+1 
Pm—\ 

Pm—n+\ 


3 . 1 . Case q-F^^\ 

(a) Setting of the Fade approximation problem 

In Table (12.21) the generating funetion and the eonstraint are established as 


(3.3) 




{aiX,a2X,asx; q)oa 
{bix,b 2 x,bjx-, q)J 


and in Table (12.51) the time evolution is ehosen as 


a\a 2 aj,q'" 

b\b 2 b-iq’^ 


(3.4) 


T : (aua 2 ,a 3 ,bub 2 ,b 3 ,m,n) (qaua 2 ,a 3 ,qbi,b 2 ,b 3 ,m,n). 


(b) Contiguity relations 

By the definition (12.91) we have the basie quantities 

(1 — bjx) 1 — byx 

(3.5) G{x) = - 4 , K{x) = - - H{x) =\ (1 - atx), 

and by the expression (12.101) we obtain the Casorati determinants 

Diix) =: cod - x/)x”’+”+iy(x)/nLd - D 2 ix) =: cix'”^”+iT(x)/(l - a,x), 

D^ix) =: Cl ^^4^(1 “ bix)il - x/g)y”+"+^y(x)/nLi(l “ 
b\ 
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where /, g, cq and ci are eonstants depending on parameters a,-, bi e C^(i = 1,2,3), m,n e Z>o 
but independent of x. Then the eontiguity relations L 2 = 0 and L 3 = 0 are expressed by 
(3.7) 

Liix) = Co(l - xf)y(x) - (1 - a 2 .^)(l - a^x)y{qx) + _ bix){l - xlg)y{x), 

bi 

L^{x) = Ci(l - xflq)y(x) + _ aix)(l - xlqg)y{x) - - b 2 x/q)(\ - b 3 xlq)y(x/q) 

bi 

where Co = co/ci and C\ = co/ci. 


Take note that in the items (c) and (d) below we study the contiguity relations L 2 = 0 and 
L 3 = 0 (13.71) for generic complex parameters = q'”, bo = q" {m, n eC^) and generic variables 
f,g (depending on parameters a,,Z 7 ; G C^,/ = 0,1,2,3) apart from the Fade approximation 
problem (12.11) with eqs. (l3.3l) and (13.41) . (see Remark 

(c) The ^-Painleve equation 

Compatibility of the contiguity relations L 2 = 0 and L 3 = 0 (13.71) gives the evolution equations 
and the constraint on the product CqCi as follows: 


(3.8) 

and 

(3.9) 


(/^-l)(/^-l) 


bob^y (/ - «2)(/ - «3)(/ - b2)(f - bs) 


aoajal (f - ai)(f - bi) 

, ^ 11 , ,(.?-! /aiiig - 1 /a3)(g - 1 /b2)(g - 1 /b^) 

(fg - 1 )(/^ - 1 ) = qaibi ---—---——--- 

{g - bi/aoa2a3){g - qbobila2a3) 


CqCi = uoibi - aoa2a3g){qbobi - a2a3g)/b\. 


The evolution equations (13.81) are equivalent to the (^r-Painleve equation of type Cg ^ given in 
iiiiibiiiiiniiii. The 8 singular points in coordinates (/, g) are on the two lines / = 00 and 
g = 00 and one curve fg=las follows: 


(/,g) = (1/(32, (32), (1/(33, a3), (1/^72, ^2), (1/^3, ^3), 

(ai, 00 ), (hi, 00 ), ( 00 , hi/(3o(32(33), ( 00 , (2'hohi/(32(33). 


(3.10) 

(d) The Lax pair 

The contiguity relations L 2 = 0 and L 3 = 0 (13.71) give two scalar Lax equations Li = 0 and 
L 2 = 0 expressed by 

(3.11) 

( 3 oho(l - bix/q)(\ - b 2 xlq){l - b^xjq) 


Li{x) = 


+ 


(1 - fx/q) 

(1 - a\x){\ - a 2 x){\ - a^x) 
q{\ - fx) 


\y{qx) 


ai(l - a 2 xlq){\ - a^xjq) 
^ 0 ^ 2 / 73(1 - bixlq){g - xlq) 
hoh 2 h 3 (l - bix){g - x) 


[y (x/q) 


}'(^)] 


-y(.^)] 


apai r/ _ h 2 h 3 g \ / _ hoh 2 ^\ 

b2b3g^\ (?<3o<3l/\ <3i / 


ai(l - a2x){l - a3x) 
xjl - a2g)(l - 6t3g)(l - b2g)il - b3g) 

(1 - fg){gq - x) 


]j(.^), 


L 2 (x) = (1 - xf)y(x) - (1 - a 2 x)(l - a 3 x)y{qx) + _ bix){\ - xlg)y{x). 

bi 
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The scalar Lax pair Li = 0 and L 2 = 0 (13.111) is expected to be equivalent to the 2 x 2 matrix 
ones in [|23l l29ll and the scalar ones in ||4l |T3l |32l by using suitable gauge transformations of 
y{x). (Note that there are some typographical errors in eqs.(30) and (31) in JH, namely the 
expressions {b 4 q' - z) and Y{x) should read {b^q' - z)t^ and Y(z) respectively.) 


(e) Special solutions 

The determinant formulae of the special solutions are given as 

1 ~ f! OLl(l ~ bj!Cl\)T[n{T^ fi+\)T (Um+Ln) 

1 ~ f! ^2 CI2 0^=l(l ~ bil(l 2 ) ^a2iT m,n+\)T( t ffi-i-X ^n) 

1 - l/(?2g _ ^2 n^=2(l “ <2i/^2) ,n+1) T })2 (t fn+ l,n) 

1 - l/bsg Z73 nL(l - ^i/h) Tg\Tm,n+l)Tb 2 irm+Un) ' 
Here the element pjt in the determinant t„, „ (13.11) is given by 


(3.12) 


(3.13) 


Pk = 




h Ik iT)( -k -^+1^3 by b2^ 

(q;q)k ^ bi b2 <33 <33 


and is the < 5 f-Appell Lauricella function (i.e. the multivariable q- hypergeometric function) 
13 defined by 


(3.14) ^X^(a,/3i,...,/3/,r;zi,...,z/) := ^ 

where |m| = mi + ... + m/. 


m ,>0 


(<T)|m|(/5l)mi • • • (fidmi jni 
(y)|m|(^)mi ■ • • (<?), 


...Z; 


mi 


Remark 7. On the transformation between the terminating and the terminating 31^2 

The element pk (13.131) expressed by the terminating (^r-Appell Lauricella series is rewritten in 
terms of the terminating ^-hypergeometric series 3^2 (12.41) (big (?-Jacobi polynomials 11^19111^) 
under the transformation of b^ = q^~^aia 2 and ai, = bib 2 . □ 

These determinant formulae of the < 2 '-hypergeometric solutions (13.121) are expected to be 
equivalent to those in dH [I3. (Note that there is a typographical error in eq.(38) in JH, 
namely the expression Ta 2 Tai{Tm,n-\) should read Ta^Ta^iTmm-i)-) Determinant formulae of q- 
hypergeometric solutions still have not been given in terms of the non-terminating ^-hypergeometric 
series 3(^2 as far as we know. 

3.2. Case q-D^^\ 


The contents of these subsections (a), (b), (c) are the same as dill- 


(a) Setting of the Fade approximation problem 


In Table (12.21) the generating function is established as 


(3.15) 


Y(x) 


(aix, a2x; q)^ 
(bix, b2x; q)oo 


and in Table (12.51) the time evolution is chosen as 


( 3 . 16 ) 


T : {ai,a2,b\,b2,m,n) {qai,a2,qbi,b2,m,n). 
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(b) Contiguity relations 

By the definition (12.91) we have the basic quantities 

2 2 
(3.17) = K{x) = \^^, H{x) = U{l-aiX), 

and by the expression (12.101) we obtain the Casorati determinants 


(3.18) 
D,{x) =: 


Co(l - x/)x”’+”+'y(x) 

n?=i(i 


D2{x) =: 


CiX 


.m+n+1 




1 - a\X 


D,(x) =: 


C2(l -(?ix)x”’+”+'y(x) 

nii(l - a,-21;) 


where /, cq, ci and C 2 are constants depending on parameters ai,a 2 , b\,b 2 6 C^, m,n e Z>o but 
independent of x. Then the contiguity relations L 2 = 0 and L 3 = 0 are expressed by 


L 2 (x) = Co(l - x/)y(x) - (1 - a 2 x)y{qx) + (1 - bix)y{x)lg, 

L 3 (x) = Ci(l - xflq)y(x) + (1 - aix)y(x)/g - - b 2 x/q)y(xlq) 


where Co = co/ci,Ci = cq/ci andg = ci/c 2 . 


Take note that in the items (c) and (d) we study the contiguity relations L 2 = 0 and L 3 = 0 
(13.191) for generic complex parameters ao = q'", bo = q'' (m, n e C^) and generic variables f,g 
(depending on parameters e C^,i = 0,1,2) apart from the Fade approximation problem 
(12.11) with eqs. (l3.15l) and (13.161) . (see Remark|^ 

(c) The ^-Painleve equation 

Compatibility of the contiguity relations L 2 = 0 and L 3 = 0 (13.191) gives the evolution equa¬ 
tions and the constraint on the product CqCi as follows: 

1 (f-ai)(f-bi) {g - bi/aoa2){g - ai/bob2) 

qaobo (/ - a2)(f - (^ 2 ) ’ ^ (g - l)(g - ^Iqaobo) 

and 

(3.21) CoC,={\-g){\-qaobog)lg\ 

The evolution equations (13.201) are equivalent to the (^r-Painleve equation of type given in 
nil Els El [13. The 8 singular points in coordinates (/, g) are on the four lines / = 0 , / = 00 , 
g = 0 and g = 00 as follows: 

(3.22) 

(/,g) = (ai,0),(Z7i,0),(0,Z7i/aoa2),(0,ai/Z7o(?2),(oo, l),(oo, l/f^aoZ^o),(« 2 , 0 °),(^ 2 , oo)- 


(d) The Lax pair 
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The contiguity relations L2 = 0 and L3 = 0 (13.191) give two scalar Lax equations Li = 0 and 
L2 = 0 expressed by 

xiai -bob2g)(bi -ao«2g)i 


= [(1 - g) (1 - qaobog) - 
g(l - aix)(l - a2x) 


(3.23) 


+ 


1 -fx 


[y(qx) - 


f 

1 - b\x 
g(l - a2x) 




T(^)] 


+ 


qaobQgjl - bixlq){\ - b2xlq) 
1 - fxlq 


\y {xjq) 


g(l -a2xlq) 
1 - b\xlq 




L 2 {x) = (1 - xf)y{x) - (1 - a 2 x)y{qx) + (1 - bix)y{x)lg. 

The scalar Lax pair Li = 0 and L2 = 0 (13.231) is equivalent to the 2x2 matrix ones in |I51 |T8]| 
and the scalar one in IfT^I^ by using suitable gauge transformations of y{x). 


(e) Special solutions 

The determinant formulae of the special solutions are given as 

\ — f!a\ FIi=l(l ~ bila\) Tai(Tm,n+l)T^^ iTm+l,n) 


(3.24) 


1 - //«2 a2 nii(l - bila2) Ta^{T, 


m,n+l 


)r '(Tm+l,n)’ 


g = 


aid - bjai) 


'^ai dm,n+\ 


)r^/(T„,+i,„) 


(7"a2(l - b2/a2) TadTm,n+\)Ta^{Tm+id 


where the element Pk Va the determinant t„, „ (13.11) is given by 


(3.25) 


Pk 


jk(^2/b2;q)k 

t>2 - 

iq', q)k 



q \ai/bi 

b2q^^^^ I a2 


bi \ 
;q,—q). 


The element (13.251) is expressed in terms of the terminating (^r-hypergeometric series 2V1 (12.41) 
(little (^r-Jacobi polynomials [0|9l|T^). These determinant formulae of the (?-hypergeometric 
solutions (13.241) are expected to be equivalent to those in [|T3]I and the terminating case of those 
in [|20l . 


3.3. Case^-Ai^\ 

(a) Setting of the Fade approximation problem 

In Table (12.21) the generating function is established as 


(3.26) 


y(x) := 


(aix, a2x; q)o 


(bix; q)oo 

and in Table (12.51) the time evolution is chosen as 

(3.27) T : {ai,a 2 ,b\,m,n) {qai,a 2 ,qbi,m,n). 


(b) Contiguity relations 
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By the definition (12.91) we have the basic quantities 

(3.28) G(x) = - - K(x) = H(x) = (1 - a,x)(l - a 2 x), 

(1 - aiA)(l - a2-^) 1 - aijc 

and by the expression (12.101) we obtain the Casorati determinants 

(3.29) 

Coil - xf)x^^'^^^Y(x) . cx^^'^^^Yix) C2A»-^"^iy(x) 

^ (1 - aiA)(l - a2x) ’ ^ l-aijc ’ ^ (1 - aix)(l - a2x) 

where /, Co,Ci and C 2 are constants depending on parameters ai,a 2 ,bi e C^, m,n e Z>o but 
independent of x. Then the contiguity relations L 2 = 0 and L 3 = 0 are expressed by 


L 2 ix) = Co(l - x/)y(A) - (1 - a 2 x)yiqx) + (1 - bix)yix)lg, 
L^ix) = Ci(l - xf/q)yix) + (1 - aix)yix)lg - q’^^"*^yixlq) 


where Q = Colc\,C\ = Cq/ci andg = Ci/c 2 . 


Take note that in the items (c) and (d) below we study the contiguity relations L 2 = 0 and 
L 3 = 0 (13.301) for generic complex parameters ao = q"\ bo = q" (m, n e C^) and generic vari¬ 
ables f,g (depending on parameters ao,a\,a 2 ,bo,b\ G C^) apart from the Fade approximation 
problem (12.11) with eqs. (13.261) and (13.271) . (see Remark!^ 


(c) The ( 2 '-Painleve equation 

Compatibility of the contiguity relations L 2 = 0 and L 3 = 0 (13.301) gives the evolution equa¬ 
tions and the constraint on the product CqCi as follows: 


(3.31) 


1 if-ai)if-bi) aia2 g-bilaoa2 

- qaobo /(/- 02 ) ’ bo (g - l)ig - 1 /quoboY 


and 

(3.32) CoCi = (1 - g)il - qaobog)lg\ 

The evolution equations (13.311) are equivalent to the ( 2 '-Painleve equation of type ^ given in 
naisiiioiiii. The 8 singular points in coordinates (/, g) are on the four lines f = 0, f = oo,g = 0 
and g = 00 as follows: 

(3.33) if,g) = iau0),ibu0),i0,bi/aoa2),ioo,l),ico,l/qaobo),ia2,c<)),ie,-ai/boe)2 

where the last point is a double point at (0, 00 ) with the gradient fg = -ailbo. (The meaning of 
the double point is also written in II 8 II.) 


(d) The Lax pair 
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The contiguity relations L 2 = 0 and L 3 = 0 (13.301) give two scalar Lax equations Li = 0 and 
L 2 = 0 expressed by 



(3.34) 


L 2 {x) = (1 - xf)y{x) - (1 - a 2 x)y{qx) + (1 - bix)y{x)lg. 


The scalar Lax pair Li = 0 and L 2 = 0 (13.341) is equivalent to the 2x2 matrix one for the 
^-Painleve equation of type q-P{AA) in IfT^ and the scalar one in lfT3]l by using a suitable gauge 
transformation of y(ji;). 

(e) Special solutions 

The determinant formulae of the special solutions are given as 



(3.35) 


q 0-2 Ta 2 iTm,n+\^T{Tm+l,n) 


where the element ptin the determinant Tm,n (13.11) is given by 



(3.36) 


The element (13.361) is expressed in terms of the terminating ( 2 '-hypergeometric series 2 V 1 (12.41) 
((^-Laguerre polynomials [ISIHIISl). These determinant formulae of the (^r-hypergeometric solu¬ 
tions (13.351) are expected to be equivalent to those in IfOll and the terminating case of those in 121. 

3.4. Case ^-(^ 2 -P 

(a) Setting of the Fade approximation problem 

In Table (12.21) the generating function is established as 

(3.37) Y(x) := (a^x, a 2 x; q)oo 
and in Table (12.51) the time evolution is chosen as 

(3.38) T : {a\,a2,m,n) {qa\,a2,m,n). 

(b) Contiguity relations 

By the definition (12.91) we have the basic quantities 



K{x) = --, H{x) = (1 - aix)(l - a 2 x), 

1 - a^x 


1 
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and by the expression (12.101) we obtain the Casorati determinants 
(3.40) 

^ - aix){\ - a2x) ’ ^ 1-aix ’ ^ (1 - aix)(l - a2-i^) 

where /, Cq, C\ and C 2 are eonstants depending on parameters ai,a 2 6 C^, m,n e Z>o but inde¬ 
pendent of X. Then the eontiguity relations L 2 = 0 and L 3 = 0 are expressed by 

L 2 (a) = Co(l - xQy(x) - (1 - a 2 x)y(qx) + y{x)/g, 

Lsix) = Ci(l - xf/q)y(x) (1 - aix)y{x)/g - q’^^"^^y{xlq) 

where Q = Co/ci,Ci = co/ciand g = Ci/c 2 . 


Take note that in the items (c) and (d) below we study the contiguity relations L 2 = 0 and 
L-i = 0 (13.411) for generic complex parameters qq = q’”, bo = q’' (m, n e C^) and generic 
variables /, g (depending on parameters ao,ai,a 2 ,bo e C^) apart from the Fade approximation 
problem (12.11) with eqs. (l3.37l) and (13.381) . (see Remark!^ 


(c) The ( 2 '-Painleve equation 


Compatibility of the contiguity relations L 2 = 0 and L 3 = 0 (13.411) gives the evolution equa¬ 
tions and the constraint on the product CqCi as follows: 


(3.42) 


gg = 


1 /-6tl 

qaobo f - a2 


// = - 


a\a2 _ g _ 

bo (g - V(g - i/^ciobo) 


and 


(3.43) CoCi = (1 - g)(l - qaobog)lg\ 

The evolution equations (13.421) are equivalent to the ( 2 '-Painleve equation of type (A 2 -i- Ai)^’\ 
namely Piv, given in |l7l[8l[l3l|2l]|. The 8 singular points in coordinates {f,g) are on the four 
lines / = 0 , / =oo,g = 0 and g = 00 as follows: 


(3.44) (/,g) = (ai,0),(oo, l),(oo, l/, 2 -ao(?o),(a 2 ,oo),( 6 ,-ai/( 7 oe) 2 ,(f,-e/< 3 o« 2 ) 2 - 


Here the fifth point is a double point at (0, 00 ) with the gradient fg = -ailbo and the sixth point 
is a double point at (0,0) with the gradient glf = -l/aoa 2 . (The meaning of the two double 
points are also written in IfSlI.l 


(d) The Lax pair 

The contiguity relations L 2 = 0 and L 3 = 0 (13.411) give two scalar Lax equations Li = 0 and 
L 2 = 0 expressed by 


Li(x) = - ^(1 - a 2 x/q)y(x)] 


(3.45) 


-I- 


-I- 

L2(x) = ( 


1 - fxjq 
g(l - aix){\ - a2x) 

1 - fx 
aoa\a2gx 


y(qx) - 


1 


/ 


-h (g - l){qaobo - 1) 


gi 

- a2x) 

1) 

y{x). 


y(^) 


xf)y(x) - (1 - a 2 x)y(qx) + y{x)/g. 































THE FADE INTERPOLATION METHOD APPLIED TO ^-PAINLEVE EQUATIONS II (DIEEERENTIAL GRID VERSIOIN^ 


The scalar Lax pair L\ = 0 and L 2 = 0 (13.451) is equivalent to the 2x2 matrix one for the 
^-Painleve equation of type q-PiA^)* in IfT^ and the scalar one in iTTSlI by using a suitable gauge 
transformation of y{x). 


(e) Special solutions 

The determinant formulae of the special solutions are given as 


(3.46) 


m+l,n) 

1 - fja2 <32 Ta2(.'^m,n+\ )r HTm+l,n)’ 


g 


Cli 7^ai(Tm,«+I (Tm+U) 

q U 2 T02^"^m,n+\yi'a2 (^»i+l,n) 

where the element Pk Va the determinant t„j „ (13.11) is given by 


(3.47) 


Pk = (-1) <?^2, 


;) (< 32 ; q)k 
(<?; q)k 


M ai X 

( 

0 


The element pk (13.471) is expressed in terms of the terminating ^-hypergeometric series 2^1 (12.41) 
(Stieltjes-Wigert polynomials [I6l|9l[l6l). These determinant formulae of the ^-hypergeometric 
solutions (13.461) are expected to be equivalent to those in [|T3ll and the terminating case of those 
inM. 


4. Conclusion 


4.1. Summary. 

In this paper for the generating function Y(x) given in Table (12.21) . we established the Fade ap¬ 
proximation problem related to the (^r-Painleve equations of type ^ and (A 2 -l- Ai)^*\ 

Then for the time evolution T given in Table (12.5L we established another Fade approxima¬ 
tion problem. By solving these problems, we derived the evolution equations, the scalar Lax 
pairs and the determinant formulae of the special solutions for the corresponding (^r-Painleve 
equations. The main results are given in Section [3l 

4.2. Problems. 

Some open problems related to the results of this paper are as follows: 

1 . In this paper by choosing certain time evolutions T, we applied the Fade method to each 

type ^ ^rid q-{A 2 + Moreover by choosing other time evolutions, 

we can perform similar computations. It may be interesting to study the relation between the 
Fade approximation method for the various time evolutions and Backlund transformations of 
the affine Weyl group, for example the case q-E^^^ of the interpolation problem in H. 

2 . By the results of this paper, it turned out that the Fade approximation method could be 

applied to the ^-Fainleve equations of type E^^\ ^ and (A 2 -l-Ai)*^^\ It may be interesting 

to study the degenerations between these results. 

3. In this paper we applied the Fade method of differential grid (i.e. Fade approximation) to 
the ^-Fainleve equations of type E^^\ A^ ^ and (A 2 -i- Ai)^^\ It may be interesting to study 
whether the Fade method can be also applied to additive difference Fainleve equations. 
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4 . It may be interesting to study whether the Fade method can be further applied to the other 
generalized Painleve systems, for example the ^-Gamier system in If2^ and the higher order 
^-Painleve system in [[24l . 
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